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We present an approximate theory of loss coefficients for modes of step- 
index fibers with various types of distortions and for fibers with lossy 
claddings. The fiber irregularities are assumed to be sinusoidal and random 
variations of the core-cladding interface. Formulas for the loss coefficients 
are presented and plotted for different values of the compound mode number 
M. For fiber lasers, we plot the loss coefficients as functions of the mirror 
tilt angles. 

We consider as an example a Nd-Y AG fiber laser with refractive index 
rii = 1.8 and a core radius of a = J+0 y.m operating at a wavelength of 
X = 1.06 n>n. For this example, we find that radiation losses are caused 
by Fourier components of fiber irregularities in the spatial wavelength 
range between 0.4 and 1.3 y.m. Intrinsic losses may be as low as 2a = 10~ 3 
cm -1 . It is thus desirable to limit scattering losses to values below 10~ 3 
cm" 1 . This requirement imposes tolerance restrictions of 0.01 fim on the 
permissible core radius fluctuations. For core radius fluctuations of this 
order of magnitude, mirror tilts should not exceed approximately 5 degrees. 
Cladding losses are not critical, but their influence on laser losses depends 
on the refractive index ratio of the core and cladding materials. Tolerable 
cladding losses may range from 10 to 300 cm~\ 

I. INTRODUCTION 

A cavity laser consists of an active medium that provides the re- 
quired gain and a (usually open) external cavity furnishing the feed- 
back for laser operation. A fiber laser also has gain and feedback, but 
instead of using the resonant modes of an open cavity it employs an 
optical fiber for guiding the radiation back and forth between the 
set of mirrors forming the cavity. 1,2 A fiber laser thus might be much 
narrower than a cavity laser since it need not allow space for the 
diffraction-limited beam to spread out in transverse direction. The 
width of the fiber laser is limited only by the loss of the fiber waveguide, 
which increases with decreasing fiber diameter. 
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In this paper we calculate the mode losses of step-index fibers and 
use them to estimate the losses of fiber lasers. The losses are caused 
by scattering from the rough fiber wall and by the presence of a lossy 
cladding. Figure 1 shows a schematic of the fiber laser. We assume 
that plane mirrors are placed at the end of the fiber that also contains 
the gain mechanism for the laser. Figure la shows a fiber laser with 
plane mirrors positioned exactly perpendicular to the fiber axis, 
whereas Fig. lb shows a laser with slightly tilted mirrors so that the 
wave inside the fiber, indicated schematically by a light ray, interacts 
more strongly with the fiber wall. For simplicity, we assume that the 
mirrors and the medium inside the fiber do not cause scattering and 
that only the fiber walls are slightly rough. We also assume that the 
fiber is surrounded by a lossy cladding that causes power loss via the 
evanescent field tail of the guided wave penetrating into the cladding. 
However, we consider these various loss mechanisms separately, one 
at a time. 

It is important to realize that wall roughness or other geometrical 
imperfections of the fiber geometry or inhomogeneities in the fiber 
material do not necessarily cause resonator losses. The electromagnetic 
field inside the cavity adjusts itself to any geometry and forms a 
normal mode. This normal mode of the cavity can be described as a 
superposition of coupled modes of the perfect waveguide. Henceforth, 
we shall refer to modes of the perfect structure as ideal modes or as 
perfect modes. The fiber imperfections provide the mechanism that 
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Fig. 1 — Schematic of the fiber laser with (a) perpendicular mirrors and (b) tilted 
mirrors. 
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couples the perfect modes together. Coupling among the guided modes 
does not introduce losses by itself. However, the perfect modes of the 
fiber suffer losses individually (in the absence of coupling). These 
losses are either caused by dissipative mechanisms in the fiber core or 
in its cladding, or they may be caused by fiber imperfections on a scale 
different from those that couple the guided modes. We may assume 
that imperfections with Fourier components of high spatial frequencies 
couple each ideal fiber mode to the radiation field outside the fiber 
and act as a loss mechanism. In addition, there will probably be im- 
perfections with large amplitudes but with low spatial frequencies 
that couple the ideal modes to each other. 

In a companion paper 3 we discuss the influence of mode coupling 
on the losses of the normal modes of the fiber cavity. We found that 
coupling among these modes increases the cavity losses compared to 
the losses of the lowest-order ideal fiber mode because neighboring 
modes with higher losses take part in the superposition field that 
forms the cavity mode. We also found that strongly coupled modes 
result in a normal mode of the cavity whose loss is an average of the 
losses of the participating coupled modes. However, not all the ideal 
modes of the fiber take part in forming the normal mode of the fiber 
cavity. Modes whose individual losses (radiation losses as well as 
dissipation losses) are relatively higher than the coupling strength to 
neighboring guided modes do not take part in the loss-averaging 
process. Since the losses of the ideal fiber modes tend to increase 
in proportion to the square of their (compound) mode number, modes 
of high order are, of necessity, much lossier than modes of low order. 
On the other hand, it is expected that the coupling strength of neigh- 
boring guided modes decreases with mode number. Consequently, 
only modes with relatively low mode numbers participate in forming 
the normal modes of the cavity. 

This theoretical expectation is confirmed by observation of laser 
radiation. 2 In fibers supporting a very large number of guided modes, 
only the modes of low order are excited as laser modes. Laser modes 
are identical with the normal modes of the resonant fiber cavity. It 
is thus clear that the loss of the normal mode of the laser (or fiber 
cavity) is an average value of the losses of the ideal modes that take 
part in forming the normal mode of the cavity. If only a very few 
fiber modes are taking part in forming the lasing mode, the loss (in 
the absence of pumping) of this cavity mode is simply the average 
loss of the few fiber modes that are effectively coupled to each other. 
In the presence of coupling among the guided modes, the loss of the 
resulting laser mode is thus somewhat higher than the loss coefficient 
of the fiber mode of lowest order, but mode coupling, even if strong 
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for the lowest-order fiber mode and its neighbor, does not increase the 
loss of the laser mode dramatically. The example studied in Ref. 3 
suggests that the loss of the laser mode may be at most an order of 
magnitude higher than the loss of the fiber mode of lowest order. We 
found, furthermore, that two modes with propagation constants /3i 
and /3 2 can couple effectively only if the relation /3 t — 02 = 2nv/L 
holds to an accuracy on the order of | on — a*| , where a indicates the 
loss coefficient, n is an integer, and L/2 is the resonator length. 

So far, we have assumed that the mirrors at the end of the fiber 
resonator are perfectly perpendicular to the fiber axis. Mirror tilt 
can be taken into account in the following way. Consider a light ray 
that propagates parallel to the fiber axis and strikes the tilted mirror 
of the resonator. After reflection, the ray impinges on the fiber wall 
at an angle that is twice the angle of the mirror tilt. Because scattering 
losses are proportional to the square of the angle between the incident 
ray and the fiber wall, it is clear that this ray, which originally traveled 
parallel to the fiber axis, suffers relatively high scattering loss. On the 
other hand, a ray that strikes the tilted mirror at normal incidence 
will strike the fiber wall at the mirror tilt angle shown in Fig. lb. 
Such a ray suffers less scattering loss. In fact, it would appear that 
the mirror tilt angle is the minimum angle at which rays passing back 
and forth through the cavity may strike the fiber wall. It is not 
obvious that there should be a ray path that closes on itself and still 
impinges at the tilted mirror at normal incidence. But the normal 
mode of a resonator has the tendency to minimize its losses. It will 
thus be composed of rays that make the lowest possible angle with 
the rough fiber walls. Consequently, we shall assume that the field 
in the resonator strikes the fiber wall at the mirror tilt angle. Instead 
of computing mode losses, we use the scattering losses of waves im- 
pinging on the rough dielectric interface at the mirror tilt angle to cal- 
culate the loss of a cavity with tilted mirrors. If both mirrors are tilted 
differently, the larger of the two angles should be used. 

We limit our discussion to fibers whose diameter is much larger 
than the wavelength of the radiation inside the fiber core. This assump- 
tion permits us to use a pseudo-plane-wave analysis. For simplicity, 
it is furthermore assumed that the refractive index difference between 
core and cladding material is so slight that reflectivity differences 
caused by polarization can be ignored; te and tm modes thus have 
the same losses. When we violate this assumption in some of our ex- 
amples, it should be remembered that our loss values apply to te 
polarization. 

Several types of wall roughness will be considered. The simplest 
imperfection is a sinusoidal variation of the fiber radius. A more com- 
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plicated wall distortion preserves the circular shape of the fiber but 
allows the diameter to vary randomly as a function of the longitudinal 
z coordinate. Finally, we consider a type of wall roughness that assumes 
that the Fourier spectrum of the wall distortion function is constant 
over all spatial frequencies of interest and that variations occur in 
both dimensions on the fiber surface with certain short correlation 
lengths. Scattering losses are expressed as functions of the amplitudes 
of the sinusoidal distortion or the variance and correlation lengths of 
the random distortion functions. Mode losses in the fiber and losses 
in the fiber cavity with tilted mirrors are considered for the case of 
scattering losses and the case of losses introduced by the lossy cladding. 

We find that cladding losses do not have a large influence on the 
wave loss in the fiber core, but scattering losses can be very serious if 
the amplitude of the wall roughness approaches the wavelength of 
the radiation. 

Spherically curved mirrors could reduce the losses of fiber lasers 
with larger diameter if they reduce the field intensity at the fiber wall. 
However, this loss reduction would work only for perfectly straight 
fibers with perpendicular mirrors and very large radii. Our estimates 
of fiber losses associated with tilted plane mirrors are equally valid 
for fiber cavities with tilted curved mirrors if the tilt angle is large 
enough. For straight fibers with perpendicular but curved mirrors, 
our loss results can be regarded as an upper limit. It should also be 
clear that mirror tilt can be translated into an abrupt tilt of the fiber 
axis. 

The analysis presented in this paper was performed to provide 
insight into the tolerance requirements of Nd-YAG fiber lasers. 2 Our 
numerical examples are thus geared to the parameters of this laser. 
The intrinsic losses of the fiber laser are on the order of 10" 3 cm" 1 so 
that additional losses caused by fiber irregularities or a lossy cladding 
should remain below this value. 

Exact loss formulas may be expressed in terms of Bessel functions 
so that their numerical evaluation becomes tedious. For this reason, 
we are here deriving simplified formulas that allow reasonable order- 
of-magnitude estimates to be readily calculated with the help of a 
simple pocket calculator. Such handy approximations are often more 
useful than the formidable exact formulas and serve the purpose of 
providing insight into the relevant variables of the problem. All our 
loss formulas are immediately applicable to optical fibers that support 
many modes. Their application to the fiber laser is straightforward if 
we can be sure that there is no additional fiber irregularity with low 
spatial frequency coupling the guided modes among each other. 
However, even if mode coupling exists, it is known from theoretical 3 
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and experimental 2 evidence that only modes of very low order par- 
ticipate in the lasing process. This information allows us to use the 
fiber loss results for the laser if we keep in mind that the loss predic- 
tion of the fiber mode of lowest order may somewhat underestimate 
the laser losses. For this reason, we base our discussion of laser losses 
on the mode with compound mode number M = 5. This loss estimate 
for the laser may, in fact, be pessimistic, but it provides the correct 
order of magnitude of the loss coefficient that may be used to derive 
tolerance requirements for the fiber laser. 

II. PLANE WAVE SCATTERING AT A PLANE INTERFACE 

We base our loss analysis on the results of plane wave scattering at 
the rough planar interface between two dielectric media, as sketched 
in Fig. 2. Our analysis uses the theory of coupled modes. In this anal- 
ysis, the incident plane wave is coupled to the continuum of modes of 
a medium that is divided into two half-spaces with a plane interface. 
The coupled mode theory is described in Ref. 4. 

To first-order perturbation theory, the scattered power is com- 
puted as follows. First, we determine the amplitudes ty(o>, a v ) of the 
continuum modes that are excited by the incident plane wave inter- 
acting with the rough interface 

/iy/2 r Lb 

dy / K jai f(y, z)dz, (1) 

-(Z„/2) JO 

where the continuous variables <r x and <r y label the continuum modes, 
Kj ai is the coupling coefficient between incident wave and continuum 
modes, and fit and p a are the propagation constants (z components of 
the propagation vectors) of incident and scattered (continuum) waves. 




INCIDENTRAY REFLECTED RAY 



Fig. 2 — Plane wave scattering at a plane, rough interface between two dielectric 
media with refractive indices n\ and ni. The y axis is directed perpendicular to the 
plane of the figure. 
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It is assumed that the rough surface with distortion function f(y, z) 
extends only over an area L y L t in the y and z directions while the 
remainder of the infinite interface is perfectly fiat. The coupling co- 
efficient is defined as 4 

K jgi - Jap (nl - «D[CrEa=o. (2) 

In this formula, a> designates the angular frequency of the light waves, 
e is the dielectric permittivity, E, is the electric vector of the incident 
and specularly reflected and transmitted waves of the perfect interface, 
while Z„ x , ay j indicates the electric field vector of the continuum mode. 
Label j designates the different types of continuum modes whose field 
expressions are given in the appendix and P in (2) is a power nor- 
malizing factor. The scattered power can now be calculated with the 
help of the formula 4 



P.c = i 3 ! f f\cj(*„ai)\*di 



*d<j y . (3) 



The summation extends over the different types of continuum modes, 
while the integration over the area S in the space a X) <r„ extends only 
over propagating continuum modes. 

With the help of the mode fields listed in the appendix, we derive 
the following expressions for the scattered power. For a sinusoidal 
corrugation of the surface 

f(y, z) = b sin flz, (4) 

we find 

p = »SJk(n\ - n\)L t L y (5) 

,c /3,(?ii sin <t>\ -f- ?i2 sm 02) 

Here b is the amplitude of the sinusoidal deflection, S z is the z com- 
ponent of the Poynting vector of the incident plane wave, and k 
and /3, are, respectively, the x and z components of the propagation 
vector of the incident plane wave in the medium with index n t whose 
magnitude is n-\k. For sinusoidal corrugation, the scattered plane 
waves are emitted in definite directions whose angles are defined by 4 



in medium 1 and by 



cob^-2*-^ (6) 

nik 



/3< - 8 /m 

cos 02 = ; — {') 



in medium 2. 

Next we list the scattering formula for scattering from a random 
corrugation. There is no variation of the surface in the y direction, 
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but the variation in the z direction with variance <r 2 is random with a 
correlation length D z that is much shorter than the wavelength of 
light. The total amount of scattered power from an area L y L z is 



with 



\7l2/ " 



P, c = - L v L z S z D z aWk {rii - ni)G\ ( ^ ) , (8) 



V(n?/n|) - 1 - | + (ni/n 2 2 ) arcsin (m/ni) 

7r[(n?/ni)- 1] 
0.7162(n 2 /m) - 0.6830 (n 8 /ni) 2 + 0.4312(n 2 /ni) 8 . (9) 



The polynomial was obtained as an empirical approximation of this 
function. Each component of the Fourier decomposition of the rough 
surface gives rise to two plane waves, one emitted into medium 1 and 
the other into medium 2. The directions of the two waves are related 
by Snell's law. If the angle (measured with respect to the surface) of 
the wave in medium 1, with the larger refractive index »i, becomes so 
small that the angle of the wave in medium 2 becomes imaginary, no 
wave can escape into medium 2; but there is still a wave emitted into 
medium 1. Equation (8) contains the large-angle contributions from 
waves emitted into both media. However, at small scattering angles 
where the wave in medium 2 disappears, the scattered wave in medium 
1 corresponds to a guided mode in a situation where medium 1 is the 
core of a fiber. Power scattered into guided mode directions is not lost, 
but becomes part of the "new" normal mode that establishes itself 
in the distorted fiber and is not included in (8). 

Finally, we list the expression for the total scattered power when 
the interface is rough in y and z dimensions. The correlation length 
(much shorter than the wavelength) of the distortion in y direction 
is D v , 

P, c = - L y L z D v D 2 a*S z KW(n\ - n&G* ( ? V (10) 

In the previous two cases, radiation was escaping only in the x, z plane. 
In the case of a truly random surface distortion, radiation escapes 
isotropically in all directions. When we apply our present results to 
the case of fiber scattering, we want to distinguish between two types 
of radiation. Any ray direction not associated with a guided mode 
belongs to either a refracting or a tunneling leaky wave. Refracting 
leaky waves leave the fiber core because they impinge on the fiber 
boundary at an angle that cannot be contained inside the fiber by 
total internal reflection. Tunneling leaky waves consist of rays that 
should be trapped inside the fiber core by means of total internal 
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reflection. 5 - 6 However, tunneling leak}' waves lose power by a mech- 
anism that causes energy to tunnel through an evanescent wave 
region outside the fiber core to an external caustic from which they 
can escape. Refracting leaky rays are very lossy and can be considered 
radiative power. Tunneling leaky rays may have very low losses in 
fibers with large core diameters and may well be part of the "new" 
normal mode of the fiber cavity. It is thus desirable to be able to dis- 
tinguish between power scattering into these two types of leaky rays. 
This distinction is made in the factor G 2 {ni/n*) appearing in (10). We 
write 

G 2 = G iT + Git. (11) 

G 2r incorporates only loss to refracting leaky rays, that is, rays scat- 
tered in those directions that, in a fiber, correspond to refracting leaky 
waves. G 2I incorporates the contribution from those scattering direc- 
tions that, in a fiber, would correspond to tunneling leaky rays. Both 
expressions could be represented in closed form but, since the closed 
form formulas would be too unwieldy, we prefer to list them in the 
form of integrals : 



(ni\ = f l "(1 ~ g) + gg>[g ~ Li ~ m/nyj dv 

/2r W J VI^T^I) (nl/n\)v + V« 2 - [1 - (nl/n!)] 



v(l - v*) + 2i> 2 V« 2 - [1 - (wj/wj)] 
(nl/n\)v + V» 2 - [1 - (n|/n?)] 
« 0.2666 (n 2 /ra x ) - 0.05359 (n 2 M) 2 + 0.3990 (n 2 M) 3 (12) 
and 

Qu = iVl-(ni/»i) - 37mi[1 -\ni/ni)l 

x f 1 Vl ~ v [ 2w + 1 - 3(ni/n?)> (13a) 

J (m/«i)t VSV — (ni/nf) 

G a tt 0.1364(n 2 /ni) + 0.7926 (m/ni)' ~ 0.2592 (wa/nO. 3 (13b) 

The polynomials are again empirical approximations. The functions 
G h G 2r , and G 2 = G 2r + G 2I are plotted in Fig. 3. As a matter of 
curiosity, we note that, disregarding the differences between functions 
G'i and 6' 2 , (8) and (10) become identical if we set the correlations 
length in the y direction equal to 

Dy = JL = *L. (14) 

Of course, this is a purely formal relationship, since (10) does not 
apply for a correlation length on the order of the wavelength. 

III. LOSS DUE TO POWER DISSIPATION IN MEDIUM 2 

In preparation for computing the fiber losses caused by a lossy jacket, 
we consider the plane wave reflection problem shown in Fig. 2 when 
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rii is lossless and the interface is perfectly plane, but medium 2 is lossy. 
The reflection coefficient for total internal reflection from the plane 
interface is expressed by the formula 7 



r = 



k + iy 



with 

and 

Using 



k — vy 
k = Vnffc 2 - /3? 
7 = V/3 2 , - n 2 2 k\ 

712 = W2r — *»« 



(15) 

(16) 

(17) 
(18) 



and the amplitude loss coefficient for plane wave propagation in 
medium 2, 

a 2 = ri2ik, (19) 



we obtain from (15), (17), and (18) 



R = |r! 2 «l - 



4<X27l2rK 



ln\-r%)W (20) 

We have assumed that uu <K n 2r and used y 2 = (n\ — nl r )W, an ap- 
proximation that holds for incident waves whose angle (with respect 
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to the interface) is far below the critical angle. When we apply our 
results to optical fibers, this assumption means that we limit ourselves 
to modes far from cutoff. The amount of dissipated power in an area 
L y L t is now 

P d = L y LJS x (l - R). (21) 

S x is the x component of the Poynting vector. We use the relation 
S x = (n/0i)S. and obtain, from (20) and (21), 

_ 2{2a i )niLyL z S z K l ,<yy. 

The real part of the refractive index has again been replaced by the 
symbol n 2 , and we used the approximation /3,- = nik. 

IV. APPLICATION TO MULTIMODE FIBERS 

The guided-mode field in optical fibers can be approximated as 4 

E y = AJ v {kt) cos v&-*: (23) 

The power density flowing in the z direction is thus given by 

& t = \—, icotf*, (24) 

2 tot 

where P is the total power carried by the guided mode. Since the 
fiber radius is a, we obtain the power density P/(wa 2 ). The factor ^ 
appearing in (24) accounts for the fact that half the total power is 
carried by a wave traveling toward the core boundary while the other 
half travels away from the boundary after reflection. The factor 
cos 2 v<(> follows in an obvious manner from (23). Averaging over the 
entire circumference of the fiber, we obtain 

S.-\— t - (25) 

The mode losses are now obtained from the plane wave formulas of 
the last two sections by identifying L y = 2va, replacing S z with 
S t of (25) and using the formula (for heat losses, P, c is replaced by Pd) 

2a - 0- (26) 

We can thus immediately compile the following list of power-loss 
coefficients for the various fiber loss mechanisms. 
Sinusoidal radius variation of amplitude b : 

_ 6V(n? - n|) m) 

2ani(wi sin 0i + ni sin 02J 

The angles <f> x and </> 2 are defined by (6) and (7). 
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Random radius variation with correlation length D z and variance <x 2 : 

2a = — Dj**k (»? - nl)G 1 ( ^ V (28) 

nia \ w 2/ 

Random surface variation with correlation length D* (formerly- 
called D y ) in the <£ direction and D z in the z direction and variance a 2 : 

2a = — DiDJ**V{n\ - n|)(7 a ( ^ V (29) 

7ra \ W2 / 

Functions Gi, (j 2 , and (7 2 r are plotted in Fig. 3. Whether G 2 or G 2r 
is to be used depends on the length and size of the fiber. If tunneling 
leaky modes are only slightly attenuated in the length of fiber under 
consideration and can be regarded as guided modes, we must use 
(?2r ; if tunneling leaky waves are very lossy, 6*2 must be used ; in inter- 
mediate cases, an average value may be appropriate. For a discussion 
of the losses of tunneling leaky waves, see Refs. 5 and 6. Finally, we 
list the power loss coefficient for a multimode fiber with lossy cladding 
(but lossless core) with cladding power-loss coefficients 2a 2 : 

2a = ?? )w * a - (30) 

nia(n\ — nlpk 3 

It remains to specify the values of k that must be inserted into formulas 
(27) through (30). In fibers supporting only one or very few guided 
modes, k would have to be obtained as the solution of the eigenvalue 
equation. However, our formulas hold only for large fibers supporting 
many modes that are mostly far from cutoff. In this case, it is possible 
to approximate k as 4 * 

* = I m - i). (3D 

The compound mode number M is a combination of the azimuthal 
mode number v and the radial mode number /x, 

M = i> + 2 M - 2,3,4, •••. (32) 

If we are interested in the losses of a fiber cavity with tilted mirrors, 
Fig. lb suggests that we use the expression 

k = n x k sin 6. (33) 

In this case, the field impinges on the fiber wall not at the natural 
mode angle applicable for perfectly straight fibers but at a larger angle 
that is imposed by the gross deformation of the fiber or mirror 



* Eq. (31) holds for small values of v. For large v, we must replace («a) 2 -* U 2 — 
[see Ref. 4, p. 90, eq. (2.5-6)] and obtain U as the solution of J V (U) = 0. 
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geometry. For fiber lasers, it seems reasonable to associate 6 with the 
mirror tilt angle. For fibers with abrupt tilts, would be the fiber tilt 
angle. 

Our derivation of formulas for the fiber loss coefficients was based 
on the properties of plane wave interaction with a plane interface. 
It is thus clear that our equations are only approximately valid. In 
particular, they do not incorporate interference effects between 
directly scattered rays and rays that leave the fiber after repeated 
reflections inside the fiber core. Such effects are particularly pro- 
nounced for scattering from purely sinusoidal core radius variations 
because, in this case, the radiation leaves at a definite angle and may 
be enhanced or reduced by interference effects. 8 Our equations give 
an average over the maxima and minima of the loss fluctuations as a 
function of scattering angle. The formulas for scattering from random 
surface effects or heat losses in the cladding are more reliable because 
diffuse scattering causes radiation to escape in all directions and inter- 
ference effects tend to cancel out and are unimportant in the case of 
power dissipation in the cladding. The formulas derived here are handy 
order-of-magnitude approximations of the precise expressions con- 
taining Bessel functions. 4,8 

V. DISCUSSION AND NUMERICAL RESULTS 

In this section, we present loss coefficients in graphic form. We 
begin with a fiber with sinusoidal core radius variations of amplitude b 
and spatial frequency ft. Scattering losses occur only if the radiation 
can escape into the cladding. The spatial frequency range that results 
in scattering losses is thus obtained from (7) as 

(ni - ni)k < ft < (ni + n 2 )k, (34) 

where we have assumed that 0,- tt wifc. If we introduce the length of 
the spatial period as A = 2ir/ft, we obtain from (34) and k = 27r/X 

x ° > A > _J^ (35) 



n\ — ni ni + n 2 

These formulas apply, of course, also to the spatial frequency range 
that contributes to random scattering, discussed below. 

Figure 4 shows curves plotted from (27). On the horizontal axis, we 
have plotted the normalized spatial frequency Q/nJc and also the 
scattering angle 4> 2 at which the radiation escapes into medium 2. 
Beyond Q/nJc = 1, the curves form the mirror image of the section 
shown in the figure and were consequently omitted. Figure 4 was 
computed for YAG with m = 1.8 and n 2 = 1. The parameter of the 
curves in Fig. 4 is the compound mode number M defined by (32). 
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Fig. 4 — Normalized scattering loss coefficient for a cavity with perpendicular 
mirrors (relative to the fiber axis) and sinusoidal core radius variation of amplitude 
6 and spatial frequency Q. M is the compound mode number. For this set of curves, 
«i = 1.8, n 2 = 1.0. 



To obtain a feeling for the magnitude of the normalized loss coefficient 
and for the tolerance requirements, we assume that the fiber reso- 
nator has an inherent loss of 2a = 10~ 3 cm -1 and a core radius of 
a = 40 /jm. Scattering loss begins to be of concern if its magnitude 
equals the already existing cavity losses. Allowing for the possibility 
that a few fiber modes of low order are tightly coupled by some fiber 
deformation of large amplitude but with a spatial frequency below 
range (34), we use an average value of 2aa 3 /b 2 = 10. If we are willing 
to tolerate a loss of 2a = 10 -3 cm -1 for a = 40 /xm, we find as the 
maximum permissible ripple amplitude the value b = 2.5 X 10 -2 /xm. 
Figure 5 shows curves that are similar to Fig. 4 except that we have 
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assumed that the YAG fiber now carries a cladding with refractive 
index ?i 2 = 1.5. It is apparent that the cladding causes a reduction of 
the scattering loss by roughly a factor of 2 so that we can now tolerate 
a ripple amplitude that is larger by v2. 

Figure 6 still describes a cavity with a fiber with sinusoidal core 
radius variation, but in this case we have assumed that the mirrors 
are tilted by an angle 0. The tilt of the mirrors causes the field inside 
the cavity to impinge on the fiber wall at an angle that is roughly equal 
to the tilt angle. It is interesting to consider the intrinsic mode angle 
to obtain a feeling for the severity of tilt angles introduced externally. 




Fig. 5— Same as Fig. 4 but with m = 1.5. 
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Fig. 6 — Normalized scattering loss coefficient for a cavity with tilted mirrors, tilt 
angle 6, and sinusoidal core radius variation of amplitude b and spatial frequency 8 
(m = 1.8 and $2 = 45 degrees). 

By equating (31) and (33), we find for the mode angle 

For a = 40 urn, X = 1.06 urn, and Wi = 1.8, we find 6 M = 0.32 
degree for the fiber mode of lowest order, with M = 2 and 6 M = 2 
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degrees for M = 10. We have assumed that mirror tilt can be con- 
trolled fairly accurately and extended our curves only to 6 = 7 
degrees. It is now* more natural to normalize the loss coefficient as 
2aa(\ /b) 2 . At a tilt angle of 6 = 5 degrees, we may expect for ni/nz 
= 1.8 the normalized loss 2aa(\o/b)°- = 0.25 according to Fig. 6. 
With 2a = 10 -3 cm -1 , a = 40 /mi and X = 1.06 /zm, we find the ripple 
amplitude b = 4.2 X 10~ 3 /mi, which is a more stringent tolerance 
requirement than the value found for straight mirrors. 

Next we consider a cavity with a fiber with randomly varying core 
radius. The case of a cavity with perpendicular mirrors is plotted from 
(28) and (9) in Fig. 7. It is assumed that the rms amplitude of the 
random core radius variation is a and that the correlation length D z 
is much shorter than the wavelength of light. If the cavity loss is an 
average value of fiber mode losses corresponding to M = 5 in Fig. 7, 
we have 2aa 3 \ /D z a 2 = 200. With a = 40 /ton and 2a = 10"" 3 cm" 1 , 
we obtain D z a 2 = 3.4 X 10 -5 /zm 3 . For want of more information, we 
assume that a = D z , so that we have D z = a = 3.2 X 10 -2 /mi. This 
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Fig. 7 — Normalized scattering loss coefficient for a cavity with perpendicular 
mirrors and random core radius variations with variance a 2 and correlation length 
D,. M is the compound mode number, rii = 1.8. 
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Fig. 8 — Normalized scattering loss coefficient for a cavity with tilted mirrors, 
tilt angle 0, and random core radius variation for n\ = 1.8. 

value is quite comparable to the value b = 2.5 pm that we found in the 
case of a sinusoidal core radius variation. 

Figure 8 gives the normalized scattering loss for a cavity with 
random core radius variation for the case of tilted mirrors. For ni = 1 
and = 5 degrees, we find from Fig. 8 2aaXg/Z>a 2 = 4. With a = 40 
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urn, we obtain a loss of 2a = 10" 3 cm" 1 for D z d 2 = 1.2 X 10" 6 Mm 8 or 
D, = a = 1.1 X 10- 2 /xm. 

Figures 9 and 10 pertain to fibers with random core-cladding interface 
perturbations. Figure 9 describes a fiber cavity with perpendicular 
mirrors. At M = 5 we find, from Fig. 9, 2aa 3 \ 2 / D *D z a 2 = 600. With 
the usual values for loss, core radius, and light wavelength, we have 
D*D Z & 2 = 1.2 X 10- 5 M m 4 or D* = D z = & = 5.9 X 10~ 2 /am. If 
we let the two correlation lengths equal the rms variation of the inter- 
face, we find that the tolerance requirements are a little less stringent 
for a totally random surface compared to a surface that maintains its 
circular cross section and only allows the radius to vary along z. 

Figure 10 shows the normalized loss coefficient for random core- 
cladding interface perturbations (in two dimensions) for a cavity 
with tilted mirrors. For a mirror tilt of 6 = 5 degrees, we obtain from 
Fig. 10 approximately 2aa)^/D i D z a 2 = 12. With a = 40 /*m, a loss 
of 2a = 10- 3 cm" 1 is obtained for D*Dr 2 = 4.2 X 10- 7 /mi 4 orD* = D t 
= a = 2.5 X 10- 2 /xm. 
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Fig. 9 — Normalized scattering loss coefficient for a cavity with perpendicular 
mirrors and random core-cladding interface perturbations with correlation length 
D+ in azimuthal direction, D, in z direction, and variance a 1 . The solid lines apply to a 
fiber whose tunneling leaky waves may be regarded as lossless guided waves; the 
dotted lines belong to the case in which tunneling leaky waves are so lossy that they 
cannot be regarded as guided waves. M is the compound mode number, n\ = 1.8. 
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The remaining figures, 11 and 12, describe the normalized loss 
coefficient for a geometrically perfect fiber core surrounded by a 
lossy cladding. The power loss coefficient of a plane wave traveling 
in the material of the cladding is 2a 2 . Figure 11 gives the mode losses 
of the fiber as a function of the normalized frequency parameter 



V = ?5 (iif- n»*. 

Ao 



(37) 




Fig. 10— Similar to Fig. 9 except that the cavity in this case has tilted mirrors with 
tilt angle 0. 
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Fig. 11 — Absorption loss coefficient of a fiber with lossy cladding (cladding loss 
coefficient 012). This set of curves applies to the case of a cavity with perpendicular 
mirrors. The normalized frequency V is defined by (37). 

For X = 1.06 nm, a = 40 ,um, m = 1.8, and n% = 1, we obtain 
V = 355 from (37). For M = 5, we obtain from Fig. 11 approxi- 
mately n&x/ntoci = 10~ 6 . The mode losses are thus much less than 
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the cladding losses. For 2a = 10 -3 cm -1 and »i/ni = 1.8, we could 
tolerate a cladding loss of 2a 2 = 1.8 X 10 3 cm -1 . If we use n% = 1.5, 
we have V = 236 leading to nia/nxx* = 4 X 10~ 6 at M = 5. With 
m/n2 = 1.2, we can now tolerate 2a 2 = 300 cm -1 . 

Figure 12 applies to a cavity with tilted mirrors and lossy jacket. 




Fig. 12 — Absorption loss coefficient of a fiber with lossy cladding (loss coefficient 
aj). These curves apply to the case of a cavity with tilted mirrors, tilt angle B. 
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Table I — Results of numerical evaluations of cladding losses. 

The cladding loss, 2a 2 , gives rise to a mode loss 

of 2a = 10- 3 cm- 1 



V 


ni/nj 


2a, 


M = 5, = 


= 5° 


355 
236 


1.8 
1.2 


1.8 X 10 3 cm~ l 
300 


56 cm -1 
10.4 



For = 5 degrees and V = 355, we find aXg(nin 2 a 2 a 2 ) = 7 X 10 -9 . 
With ni = 1.8, 712 = 1, a = 40 nm, and 2a = 10~ 3 cm -1 , we can tol- 
erate 2a 2 = 56 cm -1 . If V = 236 and ?i 2 = 1.5, we can tolerate 2a 2 
= 10.4 cm -1 . These results are summarized in Table I. 

VI. CONCLUSIONS 

We have studied the losses of fibers and fiber resonators that are 
caused by perturbations of the core-cladding interface and by absorp- 
tion losses of the cladding material. Formulas for the loss coefficient 
were derived by using plane wave techniques, and representative ex- 
amples were displayed in the form of normalized curves. The theory 
presented here is not precise, and its application to practical cases is 
hampered by lack of knowledge of coupling among the guided modes. 
We have seen in a previous paper 3 that mode coupling tends to increase 
the cavity losses above the minimum value of the fiber mode of lowest 
order. However, the loss increase due to mode coupling results only 
in an average loss of a few of the lower-order modes that are coupled 
particularly tightly. We have thus concentrated on an average loss 
corresponding to mode M = 5 when we considered explicit loss values. 
Our results are useful to gain insight into the order of magnitude of 
fiber tolerances that must be maintained and into the amount of 
cladding losses that can be tolerated. We found that the tolerances of 
core-cladding interface perturbations are on the order of 0.01 /zm, 
while cladding losses can be allowed to be as high as 10 cm -1 in the 
worst case of a cavity with mirrors tilted by 5 degrees, or as high as 
300 cm -1 in the case of a cavity with perfectly perpendicular mirrors. 
Mirror losses were lumped in with the "background losses" of the 
cavity, which were assumed to be 2a = 10~ 3 cm -1 in all the numerical 
examples we have considered. All curves independent of mirror tilt 
can be used to obtain the losses of optical fibers because they show 
plots of fiber mode losses without being tied to an application to fiber 
resonators. 
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APPENDIX 

For the calculation of scattering losses of a plane wave impinging 
on an irregular dielectric interface, we need to know the radiation 
modes of a space consisting of two dielectric media with refractive 
indices n x and n 2 separated by a plane interface. There are several 
types of such radiation modes. In each case, we list only the E z and 
H z components of their electric and magnetic fields, since all other 
field components follow from these longitudinal components by differ- 
entiation. 9 The time dependence of the modes is understood to be of 
the form 

e iat . (38) 

We place our coordinate system so that the interface between the 
media with ni and n 2 lies in the y-z plane. We assume that n x > n 2 
and let the medium with index n 2 be in the half space x > 0. 

There are radiation modes whose fields decay exponentially in 
positive x direction for x > 0. These modes can be grouped further 
into modes with E x = and H x = 0. 

(i) Evanescent modes with E x = : 

S z = Aier^e-* 1 '***™ } 

. BA . A it „ M \ for x ^ (39) 

<1)HQ<T U J 

S.-Ai (cos a x x - i j sin cj\ r*<**+i 

3C z = i ll±. Al [sin *# - - cos a x x\ e-«'*«H-M 
UflQCTy [_ (T x J 

The parameters entering these equations are related by 

nlk* = o* + el + P 2 (41) 

and 

n 2 2 k 2 = -A 2 + (rg + /3 2 . (42) 

The fields are normalized with respect to a delta function, 

i f f(S x X* v - SyKDdxdy = P8(a x - <x' x )8(<r v - tr'„), (43) 
so that we obtain for the amplitude coefficient 

2 _ 2a 2 COMQtr 2 0-gP (AAS 

Al ~ **V*B(p + *l) ' K } 
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- for x ^ 0. (40) 



with V defined by (37) and no indicating the magnetic permeability 
of vacuum. 

(ii) Evanescent modes with H x = : 

= i r ®^ A«r**(T<W>« for ^ (45) 



3C. 



aj/xo/3A 



3C Z = i 



= A 2 (cos a x x + ^ ^ sin a x x\ «r«-wrf*) 

g 4 2 (cos <r*a: - -J - sin ^ e-'<»"+w 



. n|/c 2 < 



a)/z /3A 



■ for 3 ^ 0. (46) 



Equations (41) and (42) still apply, and the amplitude coefficient is 



A\ = 



(47) 



These first two types of modes are vahd only in a limited range of 
a x and <j„ that is determined by the requirement that A, defined by 
(42), must be a positive real quantity. 

(Hi) Full standing wave modes with E x = : 

6, = A 3j (cos Px x + By sin p&)c-*<***** ) 



3C Z = % 



: fa 



WfiOVy 



A 3j (sm Px x - Rj cos p x x)e- i( - a «« +Pz ' > 



Sz = An (cos a x x + ^ Rj sin a x x\ *-*<***** 

- A 3j (sin a x x - — Rj cos a x x J «-<(*» v+fr) 



3C Z = i 



.. /3<r a 



CO/iofj 



for a: ^ (48) 



for igO. (49) 



Equation (41) applies in this case, too, but (42) is replaced by 

# = Px + crl + P. (50) 

The coefficient Rj is arbitrary, but it is convenient to choose two values 
Ri and Ri so that the two resulting modes become mutually orthogonal. 
We choose for convenience 

Rx = (51) 

and 

R 2 = oo . (52) 



The corresponding amplitude coefficients are 



A 2 — 



Z<r x a v unol 



and 



RiAm — 



2al<rlo)noP 



(53) 



(54) 
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(iv) Full standing wave modes with H x = : 
8, = A 4y (cos Px x + Si sin pjc)r*<n*+*) 
3C = -i ^^ A 4 ,(sin Pl x - Sy cos m)^*» ,w 

PxP 

• ((TyV+fc) 



for x ^ 0. (55) 



5, = A 4j (cos <r*x + — Sj sin cr^x j e - '' 

A 4 y (sin <r 2 x - ^ S; COS ^ r'""**' 

\ ™1 P* / 



= _ t .««o«r tf n? 



(T X $ 

for x ^ 0. (56) 

Equations (41) and (50) determine the relations among the compo- 
nents of the propagation constant. The ranges of a x and <j v are limited 
to the regions where p x is real and positive. This remark applies also 
to case (Hi). 

Two sets of mutually orthogonal modes result if we choose 

St - (57) 

and 

£ 2 = oo. (58) 

The amplitude coefficients are 



at. _ 2 Px oHp /eQpP /rq-, 

Ail ~ M(nlo- x + n\p x )(p + *l) K } 

and 

Q2A2 _ 27l?(T I p*VWeoffP ffi0 v 

^ 2A42 " **nlk{nl<r x + n? Pl )(/3 2 + cr*) V ' 

fo is the permittivity of vacuum. 
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